3.4.26 \(\int \frac {1}{x^4 (a+b x^4+c x^8)} \, dx\) [326]

Optimal. Leaf size=365 \[ -\frac {1}{3 a x^3}+\frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b-\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b+\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b-\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b+\sqrt {b^2-4 a c}\right )^{3/4}} \]

[Out]

-1/3/a/x^3+1/4*c^(3/4)*arctan(2^(1/4)*c^(1/4)*x/(-b-(-4*a*c+b^2)^(1/2))^(1/4))*(1-b/(-4*a*c+b^2)^(1/2))*2^(3/4
)/a/(-b-(-4*a*c+b^2)^(1/2))^(3/4)+1/4*c^(3/4)*arctanh(2^(1/4)*c^(1/4)*x/(-b-(-4*a*c+b^2)^(1/2))^(1/4))*(1-b/(-
4*a*c+b^2)^(1/2))*2^(3/4)/a/(-b-(-4*a*c+b^2)^(1/2))^(3/4)+1/4*c^(3/4)*arctan(2^(1/4)*c^(1/4)*x/(-b+(-4*a*c+b^2
)^(1/2))^(1/4))*(1+b/(-4*a*c+b^2)^(1/2))*2^(3/4)/a/(-b+(-4*a*c+b^2)^(1/2))^(3/4)+1/4*c^(3/4)*arctanh(2^(1/4)*c
^(1/4)*x/(-b+(-4*a*c+b^2)^(1/2))^(1/4))*(1+b/(-4*a*c+b^2)^(1/2))*2^(3/4)/a/(-b+(-4*a*c+b^2)^(1/2))^(3/4)

________________________________________________________________________________________

Rubi [A]
time = 0.27, antiderivative size = 365, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 5, integrand size = 18, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.278, Rules used = {1382, 1436, 218, 214, 211} \begin {gather*} \frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \text {ArcTan}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-\sqrt {b^2-4 a c}-b}}\right )}{2 \sqrt [4]{2} a \left (-\sqrt {b^2-4 a c}-b\right )^{3/4}}+\frac {c^{3/4} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \text {ArcTan}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{\sqrt {b^2-4 a c}-b}}\right )}{2 \sqrt [4]{2} a \left (\sqrt {b^2-4 a c}-b\right )^{3/4}}+\frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-\sqrt {b^2-4 a c}-b}}\right )}{2 \sqrt [4]{2} a \left (-\sqrt {b^2-4 a c}-b\right )^{3/4}}+\frac {c^{3/4} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{\sqrt {b^2-4 a c}-b}}\right )}{2 \sqrt [4]{2} a \left (\sqrt {b^2-4 a c}-b\right )^{3/4}}-\frac {1}{3 a x^3} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[1/(x^4*(a + b*x^4 + c*x^8)),x]

[Out]

-1/3*1/(a*x^3) + (c^(3/4)*(1 - b/Sqrt[b^2 - 4*a*c])*ArcTan[(2^(1/4)*c^(1/4)*x)/(-b - Sqrt[b^2 - 4*a*c])^(1/4)]
)/(2*2^(1/4)*a*(-b - Sqrt[b^2 - 4*a*c])^(3/4)) + (c^(3/4)*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTan[(2^(1/4)*c^(1/4)*x)
/(-b + Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(1/4)*a*(-b + Sqrt[b^2 - 4*a*c])^(3/4)) + (c^(3/4)*(1 - b/Sqrt[b^2 - 4*
a*c])*ArcTanh[(2^(1/4)*c^(1/4)*x)/(-b - Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(1/4)*a*(-b - Sqrt[b^2 - 4*a*c])^(3/4)
) + (c^(3/4)*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTanh[(2^(1/4)*c^(1/4)*x)/(-b + Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(1/4)
*a*(-b + Sqrt[b^2 - 4*a*c])^(3/4))

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 218

Int[((a_) + (b_.)*(x_)^4)^(-1), x_Symbol] :> With[{r = Numerator[Rt[-a/b, 2]], s = Denominator[Rt[-a/b, 2]]},
Dist[r/(2*a), Int[1/(r - s*x^2), x], x] + Dist[r/(2*a), Int[1/(r + s*x^2), x], x]] /; FreeQ[{a, b}, x] &&  !Gt
Q[a/b, 0]

Rule 1382

Int[((d_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(d*x)^(m + 1)*((a +
 b*x^n + c*x^(2*n))^(p + 1)/(a*d*(m + 1))), x] - Dist[1/(a*d^n*(m + 1)), Int[(d*x)^(m + n)*(b*(m + n*(p + 1) +
 1) + c*(m + 2*n*(p + 1) + 1)*x^n)*(a + b*x^n + c*x^(2*n))^p, x], x] /; FreeQ[{a, b, c, d, p}, x] && EqQ[n2, 2
*n] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0] && LtQ[m, -1] && IntegerQ[p]

Rule 1436

Int[((d_) + (e_.)*(x_)^(n_))/((a_) + (b_.)*(x_)^(n_) + (c_.)*(x_)^(n2_)), x_Symbol] :> With[{q = Rt[b^2 - 4*a*
c, 2]}, Dist[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^n), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), In
t[1/(b/2 + q/2 + c*x^n), x], x]] /; FreeQ[{a, b, c, d, e, n}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && NeQ
[c*d^2 - b*d*e + a*e^2, 0] && (PosQ[b^2 - 4*a*c] ||  !IGtQ[n/2, 0])

Rubi steps

\begin {align*} \int \frac {1}{x^4 \left (a+b x^4+c x^8\right )} \, dx &=-\frac {1}{3 a x^3}+\frac {\int \frac {-3 b-3 c x^4}{a+b x^4+c x^8} \, dx}{3 a}\\ &=-\frac {1}{3 a x^3}-\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^4} \, dx}{2 a}-\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^4} \, dx}{2 a}\\ &=-\frac {1}{3 a x^3}+\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b-\sqrt {b^2-4 a c}}-\sqrt {2} \sqrt {c} x^2} \, dx}{2 a \sqrt {-b-\sqrt {b^2-4 a c}}}+\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b-\sqrt {b^2-4 a c}}+\sqrt {2} \sqrt {c} x^2} \, dx}{2 a \sqrt {-b-\sqrt {b^2-4 a c}}}+\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b+\sqrt {b^2-4 a c}}-\sqrt {2} \sqrt {c} x^2} \, dx}{2 a \sqrt {-b+\sqrt {b^2-4 a c}}}+\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b+\sqrt {b^2-4 a c}}+\sqrt {2} \sqrt {c} x^2} \, dx}{2 a \sqrt {-b+\sqrt {b^2-4 a c}}}\\ &=-\frac {1}{3 a x^3}+\frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b-\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b+\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b-\sqrt {b^2-4 a c}\right )^{3/4}}+\frac {c^{3/4} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt [4]{2} a \left (-b+\sqrt {b^2-4 a c}\right )^{3/4}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [C] Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.
time = 0.03, size = 75, normalized size = 0.21 \begin {gather*} -\frac {1}{3 a x^3}-\frac {\text {RootSum}\left [a+b \text {$\#$1}^4+c \text {$\#$1}^8\&,\frac {b \log (x-\text {$\#$1})+c \log (x-\text {$\#$1}) \text {$\#$1}^4}{b \text {$\#$1}^3+2 c \text {$\#$1}^7}\&\right ]}{4 a} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[1/(x^4*(a + b*x^4 + c*x^8)),x]

[Out]

-1/3*1/(a*x^3) - RootSum[a + b*#1^4 + c*#1^8 & , (b*Log[x - #1] + c*Log[x - #1]*#1^4)/(b*#1^3 + 2*c*#1^7) & ]/
(4*a)

________________________________________________________________________________________

Maple [C] Result contains higher order function than in optimal. Order 9 vs. order 3.
time = 0.04, size = 62, normalized size = 0.17

method result size
default \(\frac {\munderset {\textit {\_R} =\RootOf \left (c \,\textit {\_Z}^{8}+\textit {\_Z}^{4} b +a \right )}{\sum }\frac {\left (-\textit {\_R}^{4} c -b \right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{7} c +\textit {\_R}^{3} b}}{4 a}-\frac {1}{3 a \,x^{3}}\) \(62\)
risch \(-\frac {1}{3 a \,x^{3}}+\frac {\left (\munderset {\textit {\_R} =\RootOf \left (\left (256 c^{4} a^{11}-256 a^{10} b^{2} c^{3}+96 a^{9} b^{4} c^{2}-16 a^{8} b^{6} c +a^{7} b^{8}\right ) \textit {\_Z}^{8}+\left (-112 b \,c^{5} a^{5}+280 b^{3} c^{4} a^{4}-231 b^{5} c^{3} a^{3}+86 b^{7} c^{2} a^{2}-15 b^{9} c a +b^{11}\right ) \textit {\_Z}^{4}+c^{7}\right )}{\sum }\textit {\_R} \ln \left (\left (\left (-1152 c^{4} a^{11}+1184 a^{10} b^{2} c^{3}-456 a^{9} b^{4} c^{2}+78 a^{8} b^{6} c -5 a^{7} b^{8}\right ) \textit {\_R}^{8}+\left (468 b \,c^{5} a^{5}-1145 b^{3} c^{4} a^{4}+933 b^{5} c^{3} a^{3}-345 b^{7} c^{2} a^{2}+60 b^{9} c a -4 b^{11}\right ) \textit {\_R}^{4}-4 c^{7}\right ) x +\left (16 a^{7} c^{5}-104 a^{6} b^{2} c^{4}+129 a^{5} b^{4} c^{3}-62 a^{4} b^{6} c^{2}+13 a^{3} b^{8} c -a^{2} b^{10}\right ) \textit {\_R}^{5}\right )\right )}{4}\) \(314\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/x^4/(c*x^8+b*x^4+a),x,method=_RETURNVERBOSE)

[Out]

1/4/a*sum((-_R^4*c-b)/(2*_R^7*c+_R^3*b)*ln(x-_R),_R=RootOf(_Z^8*c+_Z^4*b+a))-1/3/a/x^3

________________________________________________________________________________________

Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/x^4/(c*x^8+b*x^4+a),x, algorithm="maxima")

[Out]

-integrate((c*x^4 + b)/(c*x^8 + b*x^4 + a), x)/a - 1/3/(a*x^3)

________________________________________________________________________________________

Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 6664 vs. \(2 (281) = 562\).
time = 1.83, size = 6664, normalized size = 18.26 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/x^4/(c*x^8+b*x^4+a),x, algorithm="fricas")

[Out]

1/12*(12*a*x^3*sqrt(sqrt(1/2)*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 - (a^7*b^4 - 8*a^8*b^2*c +
 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^
6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)))*arc
tan(-1/2*(sqrt(1/2)*(b^14 - 16*a*b^12*c + 102*a^2*b^10*c^2 - 328*a^3*b^8*c^3 + 553*a^4*b^6*c^4 - 457*a^5*b^4*c
^5 + 152*a^6*b^2*c^6 - 16*a^7*c^7 + (a^7*b^11 - 17*a^8*b^9*c + 113*a^9*b^7*c^2 - 364*a^10*b^5*c^3 + 560*a^11*b
^3*c^4 - 320*a^12*b*c^5)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*
b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))*sqrt((b^12*c^4 - 10*a*b^10*c^5
 + 37*a^2*b^8*c^6 - 62*a^3*b^6*c^7 + 46*a^4*b^4*c^8 - 12*a^5*b^2*c^9 + a^6*c^10)*x^2 + 1/2*sqrt(1/2)*(b^18 - 1
8*a*b^16*c + 135*a^2*b^14*c^2 - 546*a^3*b^12*c^3 + 1288*a^4*b^10*c^4 - 1792*a^5*b^8*c^5 + 1421*a^6*b^6*c^6 - 5
92*a^7*b^4*c^7 + 114*a^8*b^2*c^8 - 8*a^9*c^9 + (a^7*b^15 - 19*a^8*b^13*c + 148*a^9*b^11*c^2 - 605*a^10*b^9*c^3
 + 1374*a^11*b^7*c^4 - 1672*a^12*b^5*c^5 + 928*a^13*b^3*c^6 - 128*a^14*b*c^7)*sqrt((b^12 - 10*a*b^10*c + 37*a^
2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^
2*c^2 - 64*a^17*c^3)))*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 - (a^7*b^4 - 8*a^8*b^2*c + 16*a^9
*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^6*c^6)/
(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)))*sqrt(-(b^7
 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 - (a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c +
37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^
16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)) + sqrt(1/2)*((a^7*b^17*c^2 - 22*a^8*b^15*c^3
 + 204*a^9*b^13*c^4 - 1032*a^10*b^11*c^5 + 3075*a^11*b^9*c^6 - 5417*a^12*b^7*c^7 + 5324*a^13*b^5*c^8 - 2480*a^
14*b^3*c^9 + 320*a^15*b*c^10)*x*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 -
12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)) + (b^20*c^2 - 21*a*b^18*
c^3 + 188*a^2*b^16*c^4 - 935*a^3*b^14*c^5 + 2821*a^4*b^12*c^6 - 5292*a^5*b^10*c^7 + 6083*a^6*b^8*c^8 - 4071*a^
7*b^6*c^9 + 1449*a^8*b^4*c^10 - 248*a^9*b^2*c^11 + 16*a^10*c^12)*x)*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 -
7*a^3*b*c^3 - (a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3
+ 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7
*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)))*sqrt(sqrt(1/2)*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 - (a^7
*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 -
 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c
 + 16*a^9*c^2)))/(b^12*c^7 - 10*a*b^10*c^8 + 37*a^2*b^8*c^9 - 62*a^3*b^6*c^10 + 46*a^4*b^4*c^11 - 12*a^5*b^2*c
^12 + a^6*c^13)) - 12*a*x^3*sqrt(sqrt(1/2)*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 + (a^7*b^4 -
8*a^8*b^2*c + 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5
*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c + 16*a
^9*c^2)))*arctan(1/2*(sqrt(1/2)*(b^14 - 16*a*b^12*c + 102*a^2*b^10*c^2 - 328*a^3*b^8*c^3 + 553*a^4*b^6*c^4 - 4
57*a^5*b^4*c^5 + 152*a^6*b^2*c^6 - 16*a^7*c^7 - (a^7*b^11 - 17*a^8*b^9*c + 113*a^9*b^7*c^2 - 364*a^10*b^5*c^3
+ 560*a^11*b^3*c^4 - 320*a^12*b*c^5)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c
^4 - 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))*sqrt((b^12*c^4 - 1
0*a*b^10*c^5 + 37*a^2*b^8*c^6 - 62*a^3*b^6*c^7 + 46*a^4*b^4*c^8 - 12*a^5*b^2*c^9 + a^6*c^10)*x^2 + 1/2*sqrt(1/
2)*(b^18 - 18*a*b^16*c + 135*a^2*b^14*c^2 - 546*a^3*b^12*c^3 + 1288*a^4*b^10*c^4 - 1792*a^5*b^8*c^5 + 1421*a^6
*b^6*c^6 - 592*a^7*b^4*c^7 + 114*a^8*b^2*c^8 - 8*a^9*c^9 - (a^7*b^15 - 19*a^8*b^13*c + 148*a^9*b^11*c^2 - 605*
a^10*b^9*c^3 + 1374*a^11*b^7*c^4 - 1672*a^12*b^5*c^5 + 928*a^13*b^3*c^6 - 128*a^14*b*c^7)*sqrt((b^12 - 10*a*b^
10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c
+ 48*a^16*b^2*c^2 - 64*a^17*c^3)))*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 + (a^7*b^4 - 8*a^8*b^
2*c + 16*a^9*c^2)*sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*a^3*b^6*c^3 + 46*a^4*b^4*c^4 - 12*a^5*b^2*c^5
 + a^6*c^6)/(a^14*b^6 - 12*a^15*b^4*c + 48*a^16*b^2*c^2 - 64*a^17*c^3)))/(a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2))
)*sqrt(sqrt(1/2)*sqrt(-(b^7 - 7*a*b^5*c + 14*a^2*b^3*c^2 - 7*a^3*b*c^3 + (a^7*b^4 - 8*a^8*b^2*c + 16*a^9*c^2)*
sqrt((b^12 - 10*a*b^10*c + 37*a^2*b^8*c^2 - 62*...

________________________________________________________________________________________

Sympy [F(-1)] Timed out
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/x**4/(c*x**8+b*x**4+a),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {could not integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/x^4/(c*x^8+b*x^4+a),x, algorithm="giac")

[Out]

integrate(1/((c*x^8 + b*x^4 + a)*x^4), x)

________________________________________________________________________________________

Mupad [B]
time = 5.57, size = 2500, normalized size = 6.85 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(x^4*(a + b*x^4 + c*x^8)),x)

[Out]

2*atan(-(((-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*
b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(
-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1
/4)*((x*(81920*a^15*b*c^8 + 1024*a^11*b^9*c^4 - 13312*a^12*b^7*c^5 + 62464*a^13*b^5*c^6 - 122880*a^14*b^3*c^7)
 - (-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4
 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c
 - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1/4)*(26
2144*a^17*c^8 + 4096*a^13*b^8*c^4 - 53248*a^14*b^6*c^5 + 245760*a^15*b^4*c^6 - 458752*a^16*b^2*c^7)*1i)*(-(b^1
1 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^
3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5
)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(3/4)*1i - 128*a^1
1*b*c^9 - 16*a^9*b^5*c^7 + 96*a^10*b^3*c^8)*1i + x*(8*a^10*c^10 - 4*a^9*b^2*c^9))*(-(b^11 + b^6*(-(4*a*c - b^2
)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(
1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8
 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1/4) + ((-(b^11 + b^6*(-(4*a*c - b^2)^5
)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2
) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 +
256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1/4)*((x*(81920*a^15*b*c^8 + 1024*a^11*b^9
*c^4 - 13312*a^12*b^7*c^5 + 62464*a^13*b^5*c^6 - 122880*a^14*b^3*c^7) + (-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2)
 - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*
a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^1
1*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1/4)*(262144*a^17*c^8 + 4096*a^13*b^8*c^4 - 53248
*a^14*b^6*c^5 + 245760*a^15*b^4*c^6 - 458752*a^16*b^2*c^7)*1i)*(-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^
5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c +
 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 1
6*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(3/4)*1i + 128*a^11*b*c^9 + 16*a^9*b^5*c^7 - 96*a^10*b^3*c^
8)*1i + x*(8*a^10*c^10 - 4*a^9*b^2*c^9))*(-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c
^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*
c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b
^4*c^2 - 256*a^10*b^2*c^3)))^(1/4))/(((-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2
- 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c -
 b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*
c^2 - 256*a^10*b^2*c^3)))^(1/4)*((x*(81920*a^15*b*c^8 + 1024*a^11*b^9*c^4 - 13312*a^12*b^7*c^5 + 62464*a^13*b^
5*c^6 - 122880*a^14*b^3*c^7) - (-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a
^3*b^5*c^3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5
)^(1/2) - 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 2
56*a^10*b^2*c^3)))^(1/4)*(262144*a^17*c^8 + 4096*a^13*b^8*c^4 - 53248*a^14*b^6*c^5 + 245760*a^15*b^4*c^6 - 458
752*a^16*b^2*c^7)*1i)*(-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^
3 + 280*a^4*b^3*c^4 - a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) -
 5*a*b^4*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b
^2*c^3)))^(3/4)*1i - 128*a^11*b*c^9 - 16*a^9*b^5*c^7 + 96*a^10*b^3*c^8)*1i + x*(8*a^10*c^10 - 4*a^9*b^2*c^9))*
(-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 -
a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c + 6*a^2*b^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 5*a*b^4*c*(-(4*a*c -
b^2)^5)^(1/2))/(512*(a^7*b^8 + 256*a^11*c^4 - 16*a^8*b^6*c + 96*a^9*b^4*c^2 - 256*a^10*b^2*c^3)))^(1/4)*1i - (
(-(b^11 + b^6*(-(4*a*c - b^2)^5)^(1/2) - 112*a^5*b*c^5 + 86*a^2*b^7*c^2 - 231*a^3*b^5*c^3 + 280*a^4*b^3*c^4 -
a^3*c^3*(-(4*a*c - b^2)^5)^(1/2) - 15*a*b^9*c +...

________________________________________________________________________________________